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For a difFerentiable map {x\, X2, ■ ■ ■ , Xn) — > {Xi, X2, ■ ■ ■ , Xn) that has an inverse, we show that 
there exists a Nambu-Hamiltonian flow in which one of the initial value, say Xn, of the map plays 
the role of time variable while the others remain fixed. We present various examples which exhibit 
', the map-flow correspondence. 
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In our previous paper Ij we have shown that, for a given map {x,y) — > {X,Y), there exists a Hamiltonian flow in 
which y plays the role of time variable, iff the map is differentiable and has an inverse. We would like to present 



^ ' various examples in this paper in order to clarify the correspondence. We also show that our scheme can be generalized 
' Before starting our discussion, it will be worth-while to explain the motivation of our study. Let a; be a set of 

o ; 

variables either real or complex, and f{x) a differentiable function. For a given / we can consider a sequence of the 

O 

map 

a;^/(a;)-./(2)(a,)^...^/M(a;), ;= /(/(a,)), etc.). (1) 



We call this sequence a dynamical system generated by /. 

We encounter such systems in various occasion in mathematics, physics and in the nature. All of complex dynamical 
systems are examples in mathematics. Backhmd transformations which characterize integrable systems, connection 
formulae of Stokes geometry are other examples. Renormalization group formula is an important example of dynamical 
systems in physics. 

When a set of initial values x is given, the map ^ determines values of them after m steps. From this information, 
however, it is not obvious if one can predict the values for other initial values. It will be desirable if there is a 
systematic method to analyze the initial value dependence of the map. ^ 

In order to answer this question, we have shown, in the previous paper the following proposition 
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Proposition 1 Let a differentiable and invertible map (x, y) — > {X, Y) of the plane be given and H = H{X, Y) be the 
function given by 

fx{X.Y) 

H{X,Y) = / (det J)dx, (2) 



in which 

. r d{X, Y) 
dct J = ^ ' 
cl{x,y) 

is the Jacobian of the map satisfying 

djdetJ) _^ 
dy 

then the following system of Hamilton's equations hold: 

dX dH dY dH 



dy dY' dy dX' 



Note that the initial value y of the map plays the role of time variable. The Hamiltonian is determined if the 
Jacobian of the map is given as a function of x and x is known as a function of X and Y . The dependence of [X, Y) 
on the initial value y is determined by solving the Hamilton equations (|2Jl . In the case of det J ~ 1 the area is 
preserverd through the map and the Hamiltonian turns out to be x{X^Y) itself. 

The purpose of this paper is to generalize our scheme to higher dimensional maps. We will derive naturally a 
Nambu-Hamiltonian flow with rt — 1 Hamiltonians from an n dimensional map. Another purpose of this paper is to 
present various examples to show how they fit to our scheme. 

If the initial value of one of the variables is changed the map will determine another set of final values. As far as 
the map is fixed, the change of the final values are known as well. One might think that since the map is known, 
there is no use of the Nambu-Hamilton equations. It is, however, not obvious if there exists a universal formula which 
describes the change of the final data of an arbitrary map. What we claim is that the Hamilton equation describes 
the change in the case of any two dimensional map, whereas the Nambu equation plays the role of universal formula 
in the case of higher dimensional map. 

In fact solutions to the Nambu-Hamilton equations discussed in this paper will not give any new information, but 
simply reproduce the map under certain initial conditions which must be chosen carefully. Nevertheless we believe 
that the use of our scheme will be valuable in the sense that it provides a different point of view of the map. The 
situation is similar to the relation between the Hamilton Jacobi equations and Hamilton's equations of motion. What 
we emphasize here is the correspondence between two schemes, namely, the map and the Nambu-Hamiltonian flow. 

The generalization of Proposition 1 to the higher dimensional map will be discussed in §2. Some implications of 
the map-flow correspondence will be discussed in §3. In §4 we present four examples of two and three dimensional 



maps. The correspondence between the recurrence formula and the differential equation for the Hermite polynomials 
offers an example of two dimensional map-flow correspondence. The Henon map, which is a canonical form of the 
two dimensional polynomial diffeomorphism, will be studied as an example of nonintegrable systems. The 3-point 
KdV map will provide an example of three dimensional map. Since this is a completely integrable system, we can 
reduce the number of independent dynamical variables from three to two and obtain a two dimensional flow instead 
of three dimensional Nambu flow. We will show that the two dimensional view of the flow becomes quite complicated 
compared with three dimensional one. Another example of three dimensional map is the g-difference 3-point Painleve 
IV map. The integrability of this map is guaranteed either in the continuous limit of the map or for a particular 
values of the set of parameters. 

Although the generalization to higher dimensional map is straightforward, one dimensional map must be discussed 
separately. We will show in appendix how it can be described in terms of Hamiltonian flow. 

II. NAMBU-HAMILTONIAN FLOWS 

In this section we generalize our scheme of deriving Hamiltonian systems associated with two dimensional map to 
higher dimensional ones. Since we naturally obtain the Nambu-Hamiltonian systems |^ it will be convenient to 
summarize notations briefly before we start our discussion. 

For an n dimensional map 

{XI,X2, - ■ ■ ,Xn) -> (/l,/2, ••■,/«) (4) 

we adopt, in this paper, the definition of the Nambu bracket 2] {/i, /2, ■ • • , /n} given by the Jacobian of the map 

r f f J- 1 _ d{fi, /2, • • • , /n) 

l/l, 72, ■ ■ ■ 1 Jni -TT, T- 

a{xi,X2, ■ ■ ■ ,Xn) 

Note that the following identity holds for arbitrary function g of {xi, X2, ■ ■ ■ , Xn)'. 

^{/l, /2, ■ ■ ■ , fj-l,Xl,fj + i,- ■ ■ , fnj^Z— = {/l, /2, ■ • ■ , fj-l,g, fj + 1, ■ ■ ■ , fn}- (5) 

1=1 



If {Xi, X2, • • • , Xn) are n dynamical variables, the Nambu equations of motion[^|3j are 

dXj _ d{Hi,H2,- ■ ■ ,Hn-i,Xj) . 

where Hj, j = 1, 2, • ■ • , n — 1 are Hamiltonians. The equation of motion of a function / — f{Xi, X2, • ■ ■ , Xn) follows. 



using |(SJ), as 



df AdXj df ^diHi,H2,---,Hn-i,X,) df 



rH "2^ rii Fix : ^ 2^ 



dt ^ dt dXj ^ d{Xi,X2r--,Xn) dXj 
d{Hi,H2, ■ • • , Hn-i, /) 

d{Xi,X2, • • ■ , Xn) 
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By definition of the Nambu bracket we find 

^ = 0, , = l,2,..,n-l. 

Our question is whetlier there exists a Nambu-Hamiltonian flow associated with the n dimensional map. The answer 
is given by the following proposition. 
Proposition 2 

Let a dijjerentiable and invertible map 

X = {X1,X2,- ■ ■ ,Xn) X = {Xi,X2,- ■ ■ ,Xn), (6) 

be given and Hj, j — 1,2, ■ ■ ■ ,n — 1 be functions of (Xi,X2, • • • , Xn) defined by 

j-Xr^-l(X) 

Hj=x,{X), (j = 1,2, 2), Hn^i^ (detJ)dx„_i, (7) 

in which 

, , T d{Xi,X2, • • ■ , Xn) 

det J = — 

d[xi,X2, ■■■ ,Xnj 

is the Jacobian of the map satisfying 

g(detJ) 

dXn 

Then the following system of Nambu-Hamilton's equations hold: 



dXj _ d{Hi,H2, • • • , Hn-i,Xj) 

dXn d{Xi,X2, - ■ ■ ,Xn) 



i = l,2,---,n. (8) 



Remarks 

1. The previous proposition 1 is a special case of n = 2. 

2. The initial value a;„ of the map ^ plays the role of time variable of the Nambu-Hamiltonian flow. Any other 
choice of the time variable among xj's generates a corresponding flow. 

3. By solving the set of Nambu-Hamilton's equations © we can reconstruct the map {Hi,H2, ■ ■ ■ ,i?„_i,x„) 
{Xi, X2, ■ ■ ■ , Xn) only if we choose initial conditions properly. If we consider all Nambu-Hamilton's equations 
associated with all possible choices of the 'time' variable the map is totally restored. 

Proof 

We want to have Hamiltonians -ffj's which do not depend on x„. This is true, by definition of the Hamiltonian's 
(0, in the cases of j = 1, 2, • • • , 71 — 2, i.e., 

dH-j dxi 

-r^ = y^ = 0, j = l,2,...,n-2, (9) 
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as long as xi,X2, ■ ■ ■ ,Xn~-2 are independent on x„. This is not true in the case of a;„_i, since det J, hence Hn-i, 
depends on .t„ in general. In order to eliminate the Xn dependence of Hn-i, we relate Xn-i with Xn so that the net 
dependence of Hn-i on x„ turns neutral. 

From the expression Q of the Hamiltonian i?„-i, we first notice the relations 

'■^"-1 a det J 



dHn-1 


dXn-l 


dXji 


dxji 




_ dXn-l 


dXn 


dXn 



dot ,/ 



dXn 



^ ^ , ""-1 adet J ^ 

dot J + / — dXn-l- 



Using these relations we can derive the following: 

d{Hi,H2, • ■ ■ , Hn-i 1 Xn— 1 y 



d{xi,X2, ■ ■ ■ ,Xn) 



dXn 

dHn-1 dHn-1 

dXn-1 dXn 
^ dXn-1 



det J 



dXn 

dXn-1 dHn-1 



^"-1 9 det J , 

-t: dXn-1 

OXn 

dXn-1 J ^ J dHn-1 

■ det J 



dXr^ 



Therefore dHn-i/dxn = is satisfied iff 

dXn-l 



1 d{Hi,H2,- ■ ■ .Hn-l,Xn-l) 



is true. Similarly we can derive 



dXn det J d{xi,X2T ■ ' ' i^n) 

dHn-1 dHn-1 



d{Hi,H2, • • • , Hn-i 

d{xi,X2, - ■ ■ ,Xn) 



dXn-1 dXn 
dXn ^ 

dXn-1 



det J 
det J. 



dHn-1 dxn 
dXn dXn-1 



The last equality follows to the fact dxn/dxn-i = 0. 

These results, together with the fact dxj/dxn — 0, j = 1, 2, • • ■ , n — 2, can be summaried to the formula 



dxi 



1 d{Hi,H2,---,Hn-i,Xj) 



j = l,2,---,n. 



dXn det J d{xi, X2, ■ ■ ■ , Xn) ' 

It is now straightforward to show 

dXj dxk _ 1 dXj d{Hi,H2, • • • , Hn-i,Xk) 



(10) 



dXj ^ 
dx ^ 



dxk dxn det c^Xfc d{xi, X2, ■ ■ ■ , Xn 



E 



dXj d{Hi,H2, ■ ■ ■ ,Hn-i,Xk) 



^ dxk d{Xi,X2, ■ ■ • ,X„) 

d{Hi,H2, • ■ ■ , Hn-i,Xj) 
d{Xi, X2, • • • , Xn) 



j = l,2,---,n. 
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To derive the last line we used the identity jSJ. 

III. THE MAP-FLOW CORRESPONDENCE 

Having established our proposition we would like to clarify the meaning of our results. In summary the propo- 
sition provides a view to see an arbitrary diffcrcntiable map as a dynamical flow. When {xi, X2, ■ • ■ , Xn-2) are 
given as functions of {Xi, X2, ■ • • , Xn) and the Jacobian of the map is known, we can reproduce whole of the map 
(xi, X2, • • • , Xn) {Xi, X2, • • • , Xn) by solving the Nambu-Hamilton equations under proper initial conditions. There- 
fore the Nambu-Hamiltonian flow is equivalent with the map in the sense that they are different view of one relation. 
We will present various examples of such maps in the next section. 

After m steps of the map / we obtain the map 

X={x^,X2,---,Xn) ^ X = {f["'\x)ji"'\x),---,fi"^\x)). (11) 

We can apply all of the arguments of the previous section to this map. For instance the Hamiltonians Hi, H2, • ■ ■ , Hn-2 
are given by 0, i.e., Hj = Xj = /j "''\x), j = 1, 2, • • ■ , n — 2. Only the difference from the single map is that the 
expressions of x become complicated in the case of m steps when they are written explicitly in terms of X . If j''^^ 
is the Jacobi matrix from the variable a?'^'^^ = f'^^'^{x) to cc('^+^' = f{x'^^^), the Jacobi matrix J of the map (|11|) is a 
product of J(i), J(2), . . . , Hence we have 

det J = det J(i) det J^^) • . . det J^") . 

From this information the Hamiltonian Hn-i{X) can be also calculated according to {Tj). This observation enables 
us to see the map after m steps also as a dynamical flow governed by the Nambu-Hamilton equations. Namely it 
is a flow of X such that n — 1 Hamiltonians Hi{X), H2{X), • • • , Hn-i{X) remain unchanged. It will become quite 
simple if the Jacobian is a constant. 

In the study of dynamical systems in one dimension it is well known that the nature of an orbit of a map is 
determined by local properties of the map function near the fixed and/orperiodic points of the map. There have 
been, however, known a little about the case of higher dimensional maps|j|. The main difficulty of the analysis of 
higher dimensional maps lies in its complexity. There is no systematic method, for example, to study local properties 
of the map when initial values are varied. This must be contrasted with the case of one dimensional map. Various 
notions well studied in the one dimensional case, such as chaos attractors, Julia sets etc., are not clear in the higher 
dimensional cases. 

We would like to emphasize here that the Nambu-Hamilton equations can be considered as equations which deter- 
mine local properties of the map when one of the initial values is varied. Let us rewrite the equation JSJ as 

^ . ^ n_j d{Hi,H2,- ■ ■ ,Hn-l) 
dXn d{X,,X2,---,Xj,---,Xn)' 



where Xj in this expression means that Xj is missing from (Xi, X2, ■ ■ ■ , Xn)- This formula shows that the variation 
of Xj is determined by the Jacobian of other variables. The set of equations (|12|l determine a one parameter orbit 
in the n dimensional space of {Xi, X2, ■ ■ ■ , Xn) such that Hi{X), H2(X), ■ ■ ■ , Hn-i{X) are fixed. When we study 
behaviour of the map along this orbit, other variables Xj ^ Xn are irrelevant. Moreover the set of equations l|12|) 
are sufficient to reconstruct the map if the initial conditions are chosen properly. Therefore the Nambu-Hamilton 
equations should provide a useful tool to analyse the map systematically. 

IV. EXAMPLES 

We will show, in this section, some examples of dynamical systems which exhibit the map-flow correspondence. 

A. Hermite polynomials 

Let us start from two dimensional maps. We consider first the recurrence formula for Hermite polynomials 7ik{x) 

Hk+i - xHk + kUk-i = 0, fc=l,2,3,---. (13) 

If we define new variables 

x^'^) = X, = nk/Hk-i (14) 



we obtain a map 



yik+l))- ^ik)_A_ 1 : fc = l,2,.... (15) 



The Jacobi matrix of this map is 



jik) 



(fe) 



y 



1 

1 k/iy'^^'^f 



After m steps of the map, {x,y) := (x'^'^\ y*^^^) — * {X^Y) := {x^"^\ y'^"'^^) , the Jacobian is given by 

(rn-l)! , , 

detJ:= ^ ^, 16) 

(y(m-l)y(m-2)...y(l))2 

and the Hamiltonian, according to (0), turns out to be 

m = 2: H{X,Y) ^ X{X -Yf 

m = 3: g(X,y)= (^~y^^^^^' , etc. (17) 

These data are sufficient to derive the Hamilton equations of motion (|3J| associated to the map H15|) . In the case of 
m = 2, for instance, they are 

^ ^2X{X ~Y), ^ = {Y-Xf + 2X{X~Y). (18) 
ay ay 



The equations (|18|l can be solved easily and we find 



when the initial values are chosen to satisfy 



X = X, Y = X , 

y 



X — cy^. 



Since the parameter c is arbitrary the last condition does not restrict the initial values {x,y). Thus we reproduced 
successfully the map from which we started. 

In the case of to = 3 the equations of motion are 

dX „2 4 



= -X 



dy (X - Y) 



2 ■ 



% - 4-X^ + ^^^ + 2X(y-X), (19) 



which are solved by the map 

2y 



X = x, Y ^x 



xy-l 

only if the initial values are constrained by 

c 1 

X = ^ -\ — , c : const. 
V V 

We could derive the Hamilton equations in which x plays the role of the time variable instead of y. In this case, we 
first notice the relation 

dY , ^ 
— = det J 

ox 

which will be verified by comparing the expressions in both sides. The Hamiltonian is thus given by 

H{X,Y)= detJ dy = Y + h{x). 

Here we introduced an arbitrary function h{x) of x , which does not affect the Hamilton equations derived from the 
map. (We have not used this possibility in (|17|l because it was irrelevant for the purpose of reproducing the map.) 
Using the nature of the map X = x, the Hamilton equations 

dX dH dY _ dH 
'dx ~ 'dY' lb! ~ ~dX 



associated to this system are 



dX dY_ _dhiX) 

dx " ' dx ^ dX ■ 



Solving these equations is rather trivial in this case and we obtain 

X^x Y{x) = -h{x). (21) 



The map- flow correspondence is achieved if h(x) + a; is given by the following continued fraction 

m— 1 TO— 1 TO — 2 2 1 



h{x) + X : 



y{rn,-l)(^x) X — X — —X — X 

We note that the map (|21() is nothing but the recurrence formula (|13() of the Hermite polynomials. On the other 
hand, if we substitute (fTIjl into the Hamilton equation and denote — dHm/dx, we obtain a differential- 
difference equation 

which can be also written as 

T-Cni {Ti'm+i - {fn + l)W„i) - H,n+i {Ti-'m ~ rnH„i-i) 

+m (H'jn_i - xH„i-i + H„i) Hm - m {H'^ - xHm + Wm+i) Ti-m-i = 



This is solved if 



"^m+l — ^"^m '^7717 "^m— 1 — 



TO 



or equivalently the Hermite differential equation 



x— ^-|-toH™ = (22) 



dx"^ dx 
holds. 

B. Henon map 

n . " . . 

Next we study the Henon map'5!| as another example of two dimensional map, but not integrable 6]. This is a 
canonical form of the two dimensional polynomial diffeomorphism. For arbitrary constants 6, c it is given by 



(23) 



X\_f y \ ( x \ ^( {X^-Y + c)/b 

Y ) {y^-bx + cj^ {y ) { X 

Since det J = 6, we have simply H — bx + const, which can be also written as 

H{X,Y) X^ -Y + c. (24) 
If we repeat the map we find different forms of the Hamiltonian 

TO = 3: H{XX) = {X'^ -y + ^^'^/b-Xh + cb 

TO = 4: H{X,Y)^{{X'^ -Y + cf/b^ -X + cf -{X"^ -Y + c)b + cb^, etc. 
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Now we like to know, when the Hamihonian is given by (|24|l and the Jacobian satisfies det J = 6, if the map is 
recovered. The Hamilton equations Q 

dy dy 

are solved by 

X{y)^y + a, Y{y) = y^ + 2ay + (3. 
The constants a, [3 can depend on x. Imposing the condition 

detJ = det[ J' i„ ' \= -i-ia' -I3)=b 

we obtain 

— (i — bx — c : constant. 

The Henon map 1)2311 corresponds to the case of a = and /? = —hx + c. When 6=1, the map is an area preserving 
map. 

In contrast to the previous example of Hcrmitc polynomials the Hcnon map is not integrable, hence there is no 
analytic expression for the Hamiltonian. This fact obstructs us from deriving the Hamilton equations in a systematic 
way. 

C. 3-point KdV map 

The purpose of this subsection is to illustrate the Nambu-Hamiltonian flow associated with the higher dimensional 
maps. To simplify the problem, however^^ we consider only three dimensional cases. In particular we study two 
examples, i.e., the discrete KdV equatiouj^ and the g-difference Painleve IV map. 

The KdV map 

l-^=X,+i-a;,„i, jeZ. (25) 

Xj Aj 

is known completely integrable. The general solutions are given in terms of the t function. It is a generalization of 
the Korteweg-de Vries equation to discrete time and discrete space j. To be specific we restrict our argument to three 
point KdV map [xi, X2, x^) = (x, y, z) with the periodic boundary conditions Xj — Xj+a. Then the equations become 

X X 

i-1 = Z-x, (26) 

y Y 
11 
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The integrability of this map is obvious since there are two constants of the map, which we denote by u and v : 



u — (- xyz, 

xyz 



1 1 1 

V — 1 \ h X + y + z. 

xyz 



Solving H2()|l for time evolution, we get two types of solutions^]. The map which we discuss in this paper is the one 
given by 




or conversely 



1 - 


'rxy- 


- xy'^z 


1 - 


' zx - 


- x^yz 


1 - 




- xyz"^ 


n 


V xy ^ 


' xy^z 


H 




- x'^yz 



V 1 + j/z + xyz'^ / 



Y 



V i + zx + xyz2 / 



1 - 


VXY - 


VX'^YZ 


n 


^XY - 


VX'^YZ 


1 - 


^YZ^ 


-XY'^Z 


1 - 


^YZ^ 


-XY^Z 



(27) 



(28) 



1. Three dimensional view 



We immediately find the Jacobian det J of the three dimensional map (a;,y, z) — s- {X,Y,Z) being 1, hence the 
volume is preserved under the map. According to the prescription of the previous section the Nambu Hamiltonians 
are 

Hi = X + ci, H2 = y + C2, ci , C2 : constants. 
The Nambu equations are given explicitly by 



dX 


d{Hi,H2) 


X^{l-Y^ + 2YZ + Y^Z^) 


dz 


d{Y,Z) 


{l + YZ + XY^Zf 


dY 




y2(i + 2ZX + 2XYZ'^ + 


dz 


d{Z,X) 


{l + YZ + xY'^zy 


dZ 


d{Hi,H2) 


1 + 2ZX + 2XYZ^ + 


dz 


d{X,Y) 


{l + YZ + xY'^zy 



If we take into account the fact that x,y are invariants of the Nambu flow, we see that the equations are solved by 
{X,Y,Z) of 123. 
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2. Two dimensional view 

We now like to see this dynamical system from another point of view. First we note that the map H27|l and l|28|l 
have two independent constants of map which we denote by r and s : 

r = xyz, s = (1 + xy){l + yz){l + zx). 



The previous constants u and v can be fixed if r and s are fixed. Using these constants we can reduce the number of 
independent variables. For instance, by using r we can eliminate z from the map and write H27|l as 



/ xy{l +ry + xy) \ 
r + y + rxy 

r'^ + ry + xy 
\ x{l + ry + xy) / 



(29) 



+rX + XY 
Y{l + rX + XY) 

XY{l + rX + XY) 



(30) 



\ r + X + rXY 
The result of the proposition 1 is applicable to this problem. 
The Jacobian of the transformation turns out to be 



det J : 



r + ry + xy 
x(r + y + rxy) ' 



(31) 



from which we can calculate the Hamiltonian 



H — r In a; + 7' j ln(r + ?/ + rxy) 



•In- 



X + rXY 



1 {r"^ + rX + XY){1 + rX + XY) 
~ ^ r + X + rXY 



(32) 



Y{1 + rX + XYf- r 

Contrary to the three dimensional view the variable y appears in the expression of the Hamiltonian explicitly. This 
y dependence can be eliminated if x also depends on y. That is, 

dH__dH_^ dl£dx _ ^ 
dy dy dx dy 



if 



dx xir'^ — 1)(1 + rx) 
dy rir"^ + ry + xy) 



(33) 



is satisfied. From H29|) and (|33|l we obtain 



dX 
dy 



dX ^ dXdx 
dy dx dy 



{r^x y + Ar^y^x — 2xy + 2r^xy + 2r^y^ — ry — y + 2r y + r^y + r^)x 
(r + y + rxy){r^ + ry + xy)r 

dY _ dY dY dx 
dy dy dx dy 

(1 — r^)(x^y + 2rxy + 2r^x + r^y + + r)(r + y + rxy) 
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(34) 

(35) 



x{\ + ry + xyY{r'^ + ry + xy)r 

On the other hand we can calculate ~dH/dY and dH/dX from (|32|l and see that they coincide with H34|l and H35|l . 
respectively. Thus we derived the Hamilton equations of motion Q. In this flow the Hamiltonian H of l|32|l and 
r = XY Z remain invariant. 

This example clearly shows us that the three dimensional Nambu equation is a proper way to characterize the 3d 
KdV flow. The same flow can be formulated in terms of two dimensional Hamilton equations, but we must pay the 
cost of cumbersome description. 



D. q-difference Painleve IV map 



Another example of three dimensional map, which we are going to present here, is the g-difference Painleve IV map. 
The continuous version of this map is known as the Painleve IV equation (P/y). It is completely integrable although 
there is no conserved quantity. A g-difference version of P/y has been introduced recently and discussed from the 
view point of its symmetry^]. This symmetry generates Backlund transformations and enables one to obtain some 
of special solutions. Nevertheless the integrability of the map itself is not yet obvious. 

The map we consider is given by 



^ aby 




cax 



1 - 


\-cz-\ 


- cazx 


H 


' ax - 


- abxy 


H 


- ax - 


- abxy 


1 - 


^by^ 


- bcyz 


1 - 


^by^ 


- bcyz 


1 - 


^cz-\ 


- cazx 



\ 



/ Z l + Y/b + YX/ba \ 
ca 1 + X/a + XZ/ac 

X 1 + Z/C+ ZY/cb 
abl + Y/b + YX/ba 

Y l + X/a + XZ/ac 
V Vc 1 + Z/C+ ZY/cb J 



(36) 



Here a, 5, c are some constants. If we have chosen the parameters such that a = b = c = 1 ora = fe = c = — 1, there 
exist constants of the map 

r = xyz, s — {1 + ax){l + by){l + cz) 
and hence the map becomes completely integrable. They are nothing but special cases of the 3-point Lotka-Volterra 



map discussed in j^. Otherwise the map is called g-difference Painleve IV 

This map is three dimensional. The Jacobian of the map is — (abc)'^ in this case. The Nambu Hamiltonians are 
now given by 



Hi^x + ci H2=y + C2. 
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The Nambu equations are 

dX d{Hi,H2) q{X/a){l+ X/a){l + Z/c + YZ/hc) 



dz ~ d{Y,Z) " b{l + Y/b + XY/ab){l + X/a + ZX/ca)' 

dY _ d{Hi,H2) _ q{l + Y/b){l + Z/c + YZ/bc) 

m ~ d{Z,X) ^ a(l + Y/b + XY/ab){l + X/a + ZX/ca) ' 

dZ _ d{Hi,H2) _ c^{Z/c){l + Z/c + YZ/bc) 

Iz ~ d{X, Y) ^ " (1 + X/a + ZX/ca){l + Y/b + XY/ab) ' 



which are solved by {X, Y, Z) of 



CONCLUDING REMARKS 



We have shown a correspondence between a map and a Nambu-Hamihonian flow. Examples presented in this paper 
are those of lower dimensional map. There is no difficulty to study higher dimensional maps if we do not hesitate 
lengthy-formulae. As we discussed in §3, our formulation can apply to every step of a sequence of the map. We would 
like to report some results elsewhere which can be obtained from the study in this direction. 

It is interesting that, by generalizing the map from two dimensional to higher, we have obtained the Nambu dynamics 
quite naturally. The Nambu equation, which was proposed 30 years ago'sj, was reformulated in 1993 by Takhtajanjsj 
from geometrical point of view. Quantization of the Nambu system has been attempted within the framework of 
deformation quantization^^. Recently it has attracted attention|ll( under the hope that it might provide a useful 
tool for the space-time quantization. 

The map-flow correspondence, which we discussed in this paper, is common to all differentiable maps which have 
Jacobian. Interpretation of the correspondence can be different depending on the problem. As we presented in §4, the 
Hermite polynomials can be obtained by solving either the map equations (|15|) or Hamilton equations (|19|l . They are 
equivalent to the recurrence formula (|13() and differential equations (|22|l . respectively. In this example a sequence of 
the map has a meaning of increase of the degree of polynomials. If we are considering a system of quantum oscillators, 
the map transforms an orbit to another of different quantum number, whereas the Hamilton equations determine 
orbits. This particular interpretation can be shared by all of orthogonal polynomials, since the Hermite polynomials 
are not exceptional. 

Our consideration of this paper is motivated by the observation found by one of the authors (A.S.) during the 
study of quantized Henon map to which the exact WKB method is appliedjl^. The exact WKB analysis of quantized 
Henon map is just to examine the Stokes geometry of quantum propagator which is represented by a multiple integral. 
However, since the higher dimensional exact WKB method is mainly developed for differential equations^^, so one 
has to map the problem to that of the differential equation. In order to deal with the Stokes geometry of higher 
order differential equations, it is essential to consider the bicharacteristic equation induced by the Borel transform of 
differential equation under consideration 13). A prototype of Hamiltonian formalism, which is indeed generalized in 
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this paper, appears as a bicharacteristic equation in which the initial value of the original Henon map plays a role of 
the time variable. 

Some of integrable systems have been studied within the framework of Nambu mechanics ^14^. Since we are interested 
in the correspondence between a discrete map and its associated flow we restricted our study to discrete integrable 
systems. Besides the KdV and Painleve IV maps, which we presented in this paper, there are known many discrete 
systems which are completely integrable. We have studied some of them, such as discrete Lotka-Volterra system, 
discrete time Toda lattice, discrete KP hierarchies, and found that the Jacobian of all of them is one. The map 
of such a system can be interpreted as a Backhind transformation from one equation to another which preserves 
integrability. 

We could discuss a theory of renormalization group. Renormalizing a mass m^^^ and a coupling constant g^^^ at 
a renormalization point fJ,^^\ we obtain a map (to'^', y*^^^ /i^^^) — ^ (m^^\ g^'^\ fi'^'^^). If the theory is renormalizable 
this procedure forms a group and the sequence of the map converges. Instead of repeating the map, however, there 
exists a convenient way to get useful information of the renormalized theory. Namely the scale invariance of the 
theory enables one to derive a differential equation, called the Gell-Mann Low equation, in which the initial value of 
fi becomes an independent variable. Solving this equation we can predict about phase properties of the theory. We 
would like to discuss on this topics from our point of view elsewhere. 
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Appendix : One dimensional map 

We could consider a simpler case of the map, i.e., a sequence of one dimensional differentiable map: 

qk^qk+i: fc = 0, 1,2, • • • ,m - 1 

and its inverse qk — > qk-i- This, however, does not correspond to any interesting system as far as the map is 
determined by the nearest neighbours. 

An interesting map is obtained if a change of qk is determined by the change of two neighbours. To be specific let 
us consider a map defined by 

qk = ct{qk+i) - f3{qk+2), fc = 0, 1, 2, • ■ • , m - 2 (37) 
with some differentiable functions a and /3. We assume (3 is invertible, so that the inverse map of H.S7|I exists: 

qk = P~^{a{qk-i) ~ qk-2), fc = 2, 3,---,m. (38) 
Writing them in derivative forms they are 



dk — akdk-i — bk-idk-2 
dk — o,kdk+i — Ck+idk+2 



(39) 
(40) 



where we have put dk ■= d/dqk and 



da{x) 



dx 



Jk := 



df3ix) 



x=qk 



ak := 



da{x) 



dx 



dp-\x) 



x=qk 



dx 



dx 



Ck := 



x=gk+i 



dp-\x) 



x=(3(qk+i) 



dx 



(41) 



x=(3(qk+i) 



Since g_i and qm+2 are not variables, we assume 

bo = Cm+l = 0. 

By iteration of H39() we can see a responce of qi under the variation of previous variables, say {qk-i, qk-2)'- 



where 



di = (a;a;_i 


- bi-i)di-2 - aik 


-281-3 




= Ak,idk- 


1 — bk-iAk+i,idk- 


2 






/ ak 


bk 


•■• 






1 


flfc+i bk+i 


•■• 










1 ak+2 bk+2 


••• 





Akj := det 


















••■ 


1 fli-l 


bi-1 



(42) 



V 



001 ai J 
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Similarly from H40|l we obtain 



dk = {cika.k+1 ~ Ck+i)dk+2 ~ akCk+2dk+3 
= Ak,idi+i - ci+iAk,i-idi+2 



(43) 



where 



Ak,i ■■= det 



/ cifc 1 ■■■ \ 

Cfc+i Sfe+i 1 • • • 

Cfc+2 ak+2 1 • • • 

••• ••• 





V 



ci-1 ai-i 1 
ci aj 



We note that H41|l shows ak = UkCk = Uk/bk- Hence A and A are related by 



AkA = Ck-- ■ ciA, 



k,l- 



(44) 



We want to find a responce of {qm-i,qrn) when the initial value qi is changed. We define {q,p), instead of qmiqi) 
and gm_i(gi), as dynamical variables which satisfy 



dmq{qi) = 1, d„,,+iq{qi) = 



and 



Cmdrn+lP + dm-lP = 0. 

We mention that, by the definition of the map (|37|l . qm+i is not a dynamical variable. 

We can prove the following: 
Proposition 

Let H he a function of {q, p) which is given by 



(45) 



(46) 



H 



Ci • • • Cm— id,n~ipdqo- 



Then Hamilton's equations 



hold. 
Proof 

To prove this we consider two cases 



dq _dH dp _ dH 

dqi dp ' dqi dq 



(47) 



(48) 



dl — Aim-ldm — CrnAi,m-2dm+l 

dl = Ai,mdm+1 — Cm+lAi,jn-ldm+2 — ^l.m^m+l 



(49) 
(50) 
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of 142() , and apply to q and p, respectively, to obtain 



= (51) 

dqi 

-f- = M,rndm^\V- (52) 

Similarly we apply 

dm = ^l,m(?o - feo^2,m9-l = Ax.jnQf) (53) 

dm-\ = Ai_,„_iao - 6o^2,rn-l<9_l = Ax^m-\dii (54) 

, which are obtained from (|43|) . to H and obtain 

- Ai,™ao7? (55) 
dm-iU = Ai,m-idoH. (56) 

Let us derive the first equation of H48|l from (|51|l and H5()|) . In fact the right hand side of (|56|) can be written, using 
(|?7jl . (O and ijnil) as 

rhs of IHHl) = C1C2 • • • Cm-lAim-ldm-lP 

= Ai„i-ldm-ip 

dq 



-dm-lP, 

dqi 



whereas the left hand side is 



Ihs of (O = — , 

op 

from which we obtain the first equation of l|48|l . The second equation of H48|l is derived by comparing (|52|l and l|55|l . 
Using ||17l, (gll and (O, 

rhs of = C1C2 • • ■ Cm-lAimdm~lP 

= -C1C2 • ■ • CmAimdm+lP 

= -Aund„i+ip = -rhs of lO 

hence the second equation of 148|l is proved. 

It is known that a one dimensional map equivalent with (|23|l can be also generated by the action 



■m — 1 



2 " ' ^ V3 

k=l k=l 



Ik + CQk 



The equations of motion are 

dS 



-qk+i-qk-i+2qk + qi, + c^0, /c = 1, 2, • ■ • , m - 1. (57) 
fc 
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Comparing with (|37() we find 

= 9fc + 2% + c, P{qk) = gfc, (58) 

hence 

flfe = 2gfc + 2, bfc = Cfe 1. 
If we specify a solution of (|46|l for p, so that the canonical pair is given by 

{q,P) = {qm,qm~l - Qm+l), (59) 

the Hamiltonian (|17|) turns out to be 

H = qo. 

We recall that qm+i in (|59() is not a dynamical variable but a parameter, which we denote as qm+i = a- The equation 
(|37|l enables us to express H in terms of . For instance 

m = 2: i/(g,p) = (p + a + 1)2 - g - 1, 
m = 3: i/(g,p) = ((p + a+l)2-g)2-i, 

where we have put c — 0. 

The correspondence between the two sequences (|23f) and (|58f) of the Hcnon map is established if we relate the 
variables by 

{X,Y) = {q + l,p + a + l). 
The generating function of this transformation is 

W{X,p) = {p + a + l){X-l). 



Before closing, we like to mention a small change of our choice of canonical variables from those of ref.[l2|- In the 
latter the canonical pair of variables is chosen as 

{(l',P') = {-qrn,qm- qrn-l) = {-q,q-p~a), 

while the Hamiltonian is given by 

dqi 

These two variables are related by the canonical transformation generated by the following function: 



W{q,p',qi) = (^^Q-P - 9 - qoqi- 



